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We report on recent progress in testing the factorization formalism of nonrelativistic 
quantum chromodynamics (NRQCD) at next-to- leading order (NLO) for J/i/> yield and 
polarization. We demonstrate that it is possible to unambiguously determine the leading 
color-octet long-distance matrix elements (LDMEs) in compliance with the velocity scal- 
ing rules through a global fit to experimental data of unpolarized J/i/i production in pp, 
pp, ep, 77, and e+e~ collisions. Three data sets not included in the fit, from hadropro- 
duction and from photoproduction in the fixed-target and colliding-beam modes, are 
nicely reproduced. The polarization observables measured in different frames at DESY 
HERA and CERN LHC reasonably agree with NLO NRQCD predictions obtained using 
the LDMEs extracted from the global fit, while measurements at the FNAL Tevatron 
exhibit severe disagreement. We demonstrate that alternative LDME sets recently ob- 
tained, with different philosophies, in two other NLO NRQCD analyses of J/ip yield and 
polarization also fail to reconcile the Tevatron polarization data with the other available 
world data. 
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1. Introduction 

While the overly successful experiments at the LHC are exploring the Higgs sec- 
tor and are systematically searching for signals of physics beyond the standard 
model (SM), we must not be carried away losing track of a longstanding, unre- 
solved puzzle in quantum chromodynamics (QCD), the otherwise well-established 
SU(3) gauge theory of the strong interactions, right in the core of the SM. In fact, 
despite concerted experimental and theoretical efforts ever since the discovery of the 
J/ip meson in the November revolution of 1974 (The Nobel Prize in Physics 1976), 
the genuine mechanism underlying the production and decay of heavy quarkonia. 
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which are QCD bound states of a heavy quark Q — c,b and its antiparticle Q, has 
remained mysterious. 

Nonrelativistic QCD (NRQCDPendowed with an appropriate factorization the- 
orem, which was conjectured in a seminal work by Bodwin, Braaten, and LepagJ^l 
and exphcitly proven through next-to-next-to-leading order for large transverse mo- 
menta p^p^ arguably constitutes the most probable candidate theory at the present 
time. This implies a separation of process-dependent short-distance coefficients, to 
be calculated perturbatively as expansions in the strong-coupling constant as , from 
supposedly universal long-distance matrix elements (LDMEs) , to be extracted from 
experiment. The relative importance of the latter can be estimated by means of 
velocity scaling rules,!^ which predict each of the LDMEs to scale with a definite 
power of the heavy-quark velocity v in the limit u <C 1. In this way, the theoretical 
predictions are organized as double expansions in as and v. A crucial feature of 
this formalism is that the QQ pair can at short distances be produced in any Fock 
state n = 2-5-1-1 ^fH]^ definite spin S, orbital angular momentum L, total angular 
momentum J, and color multiplicity a = 1, 8. In particular, this formalism predicts 
the existence of intermediate color-octet (CO) states in nature, which subsequently 
evolve into physical, color-singlet (CS) quarkonia by the nonperturbative emission 
of soft gluons. In the limit v ^ 0, the traditional CS model (CSM) is recovered in 
the case of S'-wave quarkonia. In the case of J/ijj production, the CSM prediction is 
based just on the ^s\^^ CS state, while the leading relativistic corrections, of relative 
order ©(f"*), are buih up by the ^S^o\ ^sf\ and (J = 0, 1,2) CO states. 

The CSM is not a complete theory, as may be understood by noticing that 
the next-to-leading-order (NLO) treatment of P-wave quarkonia is plagued by un- 
canceled infrared singularities, which are, however, properly removed in NRQCD. 
This conceptual problem cannot be cured from within the CSM, neither by pro- 
ceeding to higher orders nor by invoking kx factorization etc. In a way, NRQCD 
factorization}-^ appr opriate ly improved at large values oipx by systematic expansion 
m powers of m|/p|)SEIH] is the only game in town, which makes its experimental 

verification such a matter of paramount importance and general interest.^ 

The experimental test of NRQCD factorizatiorP has been among the most ur- 
gent tasks on the agenda of the international quarkonium conimunitjEl for almost 
two decades and, with high-quality data being so copiously harvested at the LHC, 
is now more tantalizing than ever. In the following, we discuss the present status of 
testing NRQCD factorization in charnionium production. 



2. Global fit to measurements of unpolarized J/if> yields 

We consider the inclusive production of J/ip mesons in collisions of two particles 
A and B. Owing to the factorization theorems of the QCD parton model and 
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Table 1. NLO NRQCD fit results for the J/ip CO LDMEs.CSl Subtracting 
from the data the estimated contributions from the feed-down of heavier char- 
monia, which are not included in the calculations, improves the quality of the 
fit. 





set A: unsubtractcd 


set B; subtracted 


^d.o.f. 


(4.97 ± 0.44) X 10-2 GeV^ 
(2.24 ±0.59) X 10-3 GeV^ 
(-1.61 ±0.20) X 10-2 GeV^ 
4.42 


(3.04 ± 0.35) X 10-2 GeV^ 
(1.68 ± 0.46) X 10-3 GeV^ 
(-9.08 ± 1.61) X 10-3 GeV^ 
3.74 



NRQCD,!2!the cross section is calculated as 
da{AB J/'ip + X) = ^ J dxidx2dyidy2 fi/A{xi)fk/i{yi)f]/B{x2)fi/j{y2) 

i.j,k,l.n 

X {0^''^[n])d(j{kl -> cc[n] + X), (1) 

where fi/A{xi) is the parton distribution function (PDF) of parton i = g,<l,q 
in hadron A — p,p or the flux function of photon i = 7 in charged lepton 
A — , e+, fk/i{yi) is SikS{l — yi) or the PDF of parton k in the resolved photon i, 
daikl — > cc[rt] + X) arc the partonic cross sections, and (O"^/''' [n]) are the LDMEs. 
In the fixed-flavor-nuniber scheme, we have q = u,d, s. In the case of e~^e~ annihi- 
lation, all distribution functions in Eq. ([Ij are delta functions. The hadronic system 
X always contains one hard parton at leading order (LO) and is t aken t o be void 
of heavy flavors, which may be tagged and vetoed experimentally.ll^ilU xhe par- 
tonic cross sections appropriate for the direct production of unpolarized J ftp mesons 
were calculated at NLO in NRQCD in Refs. [121 HSI for direct photoproduction, in 
Refs. [Ml [151 [IB [m for hadroproduction, and in Ref. [18] for resolved photoproduc- 
tion, two-photon scattering involving both direct and resolved photons, and e^e~ 
annihilation. 

In our numerical analysis, we set rric = 1.5 GeV, adopt the values of me, a, and 
the branching ratios B{J/ijj — > e+e^) and B{J/ip — > ^1^ ji^) from Rcf. 19, and use 
the one- loop (two- loop) formula for a'i^^\p), with Uf — A active quark flavors, at LO 
(NLO). As for the proton PDFs, we use set CTEQ6L1 (CTEQ6M)20 at LO (NLO), 
which comes with an asymptotic scale parameter of Aq^j-, — 215 MeV (326 MeV). 
As for the photon PDFs, we employ the best-fit set AFG04_BF of Ref. [211 We 
evaluate the photon flux function using Eq. (5) of Ref. [121 with the upper cutoff on 
the photon virtuality chosen as in the considered data set. As for the CS LDME, 
we adopt the value (©"'/'''(^S'W)) = 1.32 GeV'"^ from Ref. [231 Our default choices 
for the renormalization, factorization, and NRQCD scales are Hr = IJ-f = and 
Ha = Wc, respectively, where rriT = \/pfr+4mf is the J/ip transverse mass. The 
bulk of the theoretical uncertainty is due to the lack of knowledge of corrections 
beyond NLO, which are estimated by varying fir, jjtf, and /xa by a factor 2 up and 
down relative to their default values. 
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In Ref. 1181 we performed a global fit to high-quality data of inclusive unpolarized 
J/tp production, comprising a total of 194 data points from 26 data sets. Specifically, 
these include d pt di stributions in hadroproduction from PHENIX^^ at RHIC, CDF 
at Tevatron i^SMI and 11,123 ATLAS ,-22^ CMS ,20! ALICE, and LHCb^^l at 
the LHC; p"^, W, and z distributions in photoproduction from hJ^SI and ZEU^^H 
at HERA I and hJ^^ at HERA II; a p^ distribution in two-photon scattering 
from DELPnf^ at LEP II; and a total cross section in e^e annihilation from 
Bell^^ at KEKB. Denoting the photon, proton, and J/ip four-momenta by p^, 
Pp, and respectively, W — ^ {p^ + ppY is the 7p center-of-mass energy and 

■2 = {pj/ip ■ Pp)/{p-t ■ Pp) is the inelasticity variable measuring the fraction of the 
photon energy passed on to the J/?/" meson in the proton rest frame. We excluded 
from our fit all data points of two-photon scattering with pT < 1 GeV and of 
hadroproduction with pT < 3 GeV, which cannot be successfully described by our 
fixed-order calculations as expected. The fit results for the CO LDMEs obtained at 
NLO in NRQCD with default scale choices are collected in Table [TJ They depend 
only feebly on the precise locations of the pt cuts. 

Our calculations refer to direct J/ijj produ ction, as the da ta from Tevatron j^5 | 26 | 
do, while the data from Tevatron 11,^11 LHC,l2 g | 29 | 30 | 31 | 3 2 | 

and KEKI^lQl comprise 

prompt events and those from RHIC,'^ 

HERA,|33I31I35| 

and LEP II even non- 
prompt ones. The fit results obtained neglecting the effects due to these admixtures 
are listed in the second column of Table [T] (set A). However, the resulting error is 
small against our theoretical uncertainties and has no effect on our conclusions. In 
fact, the fraction of J/Tp events originating from the feed-down of heavier charmonia 
only amounts to about 36% for hadroproductionp^ 15% for photoproduction at 
HERA,'^9% for two-photon scattering at LEP II,E3and 26% for e+e~ annihilation 
at KEKB,"^^ and the fraction of J/^p events from B decays is negligible RHIC, 
HERA,l25l 

and LEP 11^ energies. Refitting the data with the estimated feed-down 
contributions subtracted yields the values listed in the third column of Table [1] (set 
B). The values per data point achieved by the two fits, which are specified as 

o f ii^ Table [U are to be taken with a grain of salt, since they do not take into 
account the theoretical uncert aint ies, which exceed most of the experimental errors. 

The fact that the global filfiS! successfully pins down the three CO LDMEs as 
it does is quite nontrivial by itself and establishes their universality, the more so 
as the long-standing difficulty of NRQCD to describe the photoproduction data at 
large values of z is overcome. Furthermore, their values are of order ©(w*) with 
respect to the CS LDME {0'^^'^{^s[^^))^m comphance with the velocity scaling 
rules.l^ Both observations consolidate the validity of NRQCD factorization as far as 
the unpolarized J/ip yield is concerned. 

In Fig. [U all data sets fitted to are compared with our default NLO NRQCD 
results (solid lines). For comparison, also the default results at LO (dashed lines) 
as well as those of the CSM at NLO (dot-dashed lines) and LO (dotted lines) are 
shown. In order to visualize the size of the NLO corrections to the hard-scattering 
cross sections, the LO predictions are evaluated with the same LDMEs and PDFs. 
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The yellow and blue (shaded) bands indicate the theoretical errors on the NLO 
NRQCD and CSM results. We observe from Fig. [T] that the experimental data are 
nicely described by NLO NRQCD, being almost exclusively contained within its er- 
ror bands, while they overshoot the NLO CSM predictions typically by 1-2 orders of 
magnitude for hadroproduction and a factor of 3-5 for photoproduction. In contrast 
to the LO analysis of Ref. [39l the DELPHI datcP^tend to systematically overshoot 
the NLO NRQCD result, albeit the deviation is by no means significant in view of 
the sizable experimental errors. This may be attributed to the destructive interfer- 
ence of the i^^l and contributions, which is a genuine NLO phenomenon. We 
have to bear in mind, however, that the DELPHI measurement comprises only 16 
events with > 1 GeV and has not been confirmed by any of the other three LEP II 
experiments. The Belle measurement, a{e~^e~ — J/tp -\- X) = (0.43 ± 0.13) pb,^^ 
is compatible both with the NLO NRQCD and CSM resuhs, (O.TO:^" :^^) pb and 
(0.24 ^g gg) pb , r cs p c c 1 1 vcly . However, the measured cross section was actually ob- 
tained from a data sample with the multiplicity of charged tracks in the events 
being larger than four, and corrections for the effect of this requirement were not 
performed, so that the value quoted in Ref. [10]just gives a lower bound on the cross 
section. 



3. Further tests of NRQCD factorization in unpolarized J/ip 
production 

Three data sets not included in the global fit ,^1^ from hadroproduction and from pho- 
toproduction in the fixed-target and colliding-beam modes, are nicely reproduced 
by our NLO NRQCD predictions, as may be seen from Figs. [2] and |3l They were 
taken by the ATLAS Collaboratioil^ni at the LHC, by Denby et al.SB at the Fermi- 
lab Tagged-Photon Spectrometer, and by the ZEUS CoUaboratioiI^l at HERA II. 
The Xd o f values evaluated using our default NLO NRQCD predictions read 10.740 
0.40, and 7.50, respectively. We conclude that NRQCD factorization passes this non- 
trivial test, which, in the case of Refs. HOI SIl probes kinematic regions far outside 
those covered by the global fit.^lSl 



4. J/'0 polarization 

The polarization of the J/^ meson is conveniently analyzed experimentally by 
measuring the angular distribution of its leptonic decays, which is c usto marily 
parametrized using the three polarization observables Ag, A^, and Ae^, aJ^Sl 

W{0, (f)) xl + Xg cos^ 9 + X4, sin^ e cos(20) -t- Ag^ sin(26') cos (2) 

where 6 and (f> are respectively the polar the azimuthal angles of Z"*" in the J/f/; 
rest frame. This def inition depend s on the choice of coordinate frame. In the ex- 
perimental analyses j '^^l^^ l ^^l^^l^'^ l the helicity (recoil), CoUins-Soper, and target 



^This value is reduced to 4.88 if the data point at the largest value of px is omitted. 
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frames were employed, in which the polar axes point in the directions of ~{pp+Pp), 
Pp/lPpl ~ Pp/\Pp\j —pp, respectively. The values Xe = 0,+l,— 1 correspond 
to unpolarized, fully transversely polarized, and fully longitudinally polarized J/^/) 
mesons, respectively. The alternative notation A = Ag, /x = Ag^, and v = 2X^ is 
frequently encountered in the literature. In Refs. |45l |46l Xg is called a. 

Working in the spin density matrix formalism and denoting the z component of 
S by i, j = 0, ±1, we have 

daii-daoo . dai,-i V2Redaio 

dan + daoo dan + da^o dan + daoo 

where daij is the ij component of the differential cross section. An expression of 
daij in terms of PDFs and partonic spin density matrix elements may be found in 
Eq. (3) of Ref. M 

Our results for direct photoproductiorPS! are shown in Fig. 31 We compare our 
NLO predictions for the parameters A and v as functions of px and z with mea- 
surements by the HI CoUaboratiorPSl in the helicity and CoUins-Soper frames and 
by the ZEUS CollaboratioiPS in the target frame. Unfortunately, the IliSS and 
ZEU^^ll data do not yet allow us to distinguish the production mechanisms clearly. 
However, kinematical regions can be identified in which a clear distinction could 
be possible in more precise experiments at a future ep collider, such as the CERN 
LHeC.49 At higher values of pT, NRQCD predicts the J/-0 meson to be largely un- 
polarized, in contrast to the CSM. In the z distributions, the scale uncertainties are 
sizable, and the error bands of the CSM a nd NR QCD predictions largely overlap. 

Our results for direct hadroproductiorlSSISlI are shown in Fig. [5l We compare 
our predictions for the parameters Xe and A^ as funct ions of pt in the helicity and 
Collins-Soper frames with the measurements by CDpSSESI and ALICE.E3 In the 
helicity frame, the CSM predicts the J/ip polarization to be strongly longitudinal 
at NLO, while NRQCD predicts it to be strongly transverse. In the Collins-So per 
frame, the situation is inverted. The precise CDF measure men t at Tevatron run 11,^^ 
which is partially in disagreement with the one at run I,S51 found the J/'tp mesons 
to be largely unpolarized in the helicity frame, which is in contradiction with both 
the CSM and NRQCD predictions at NLO. The early ALICE datcPS is, however, 
compatible with NRQCD at NLO, favoring NRQCD over the CSM. 



5. Comparisons with the literature 

After our NLO NRQCD studies of J ftp polarizationj ^^ l ^Q ^^ two others appeared, 
which are, however, limited to hadroproduction. In Ref. I52| it was shown that the 
measured hadroproduction cross sections and the CDF II polarization measure- 
ment can be simultaneously described by NRQCD at NLO with one of the three 
CO LDME sets listed in the fourth column of Table [2j In Ref. |53l the polarization 
of promptly produced J/V' mesons was studied by also including the feed-down 
from polarized XcJ and ip' mesons as described in Refs. [Ml |55l To this end, the 
CO LDMEs of the Xcj and ip' mesons were fitted to LHCb (and CDF) unpolarized 
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production data, and the resulting cascade decay rates into J/ mesons were then 
used as feed-down contributions to determine the J/ip CO LDMEs in a fit to unpo- 
larized J/tp production data from LHCb and CDF with pr > 7 GeV. The resulting 
LDMEs may be found in the third column of Table [2] Reference [53] predicts the 
J/V' polarization to be moderately transverse in the helicity frame. 

In Fig. ini we systematically compare the analyses of Refs. |50l [511 [52l [53|as rep- 
resented by the CO LDME sets in Table [2| with regard to their performances in 
describing the unpolarizedJ7'0 yields measured in e~^e~ annihilation by Belle in 
photoproduction by Hlj ^^ l ^^ land in hadroproduction by CDF li23and ATLAS,S2l 
as well as the J/^ polarization observable Xe in the helicity frame as measured 
by CDF II.SSl We observe that none of the LDME sets can describe all the data 
sets. While the CO LDMEs of Ref. [T8| yield a good description of the unpolar- 
ized J/ip yields, there is a strong disagreement with the CDF II measurement of 
J/ip polarization. On the other hand, the CO LDMEs of Ref. [52| can describe all 
hadroproduction data, but lead to overshoots by factors of 4-6 for e+e" annihi- 
lation and photoproduction. Finally, the CO LDMEs of Ref. [53| yield predictions 
which, in all cases, fall between those of the other two options. 



Table 2. LDME sets determined in Refs. [T8ll52l [53l and used in Fig. |6] In Ref. [52] two alter- 
native sets are provided besides the default one. The analyses of Refs. l50l 1511 [52l onlv refer to 
direct J/xp production. 





Butenschoen, 


Gong, Wang, 


Chao, Ma, Shao, Wang, 


ZhanJSH 




KniehlSl 


Wan, Zhan#Sl 


default set 


set 2 


set 3 




1.32 GeV^ 


1.16 GeV^* 


1.16 GcV« 


1.16 GcV^ 


1.16 GeV^* 




0.0497 GeV=^ 


0.097 GeV^* 


0.089 GeV^ 





0.11 GeV^ 




0.0022 GeV^ 


-0.0046 GeV^ 


0.0030 GeV^ 


0.014 GeV^ 







-0.0161 GeV^ 


-0.0214 GeV^ 


0.0126 GeV^ 


0.054 GeV^ 





(0''''(3sf' )) 




0.758 GeV^ 












-0.0001 GeV^ 








^0V'(35l81)^ 




0.0034 GeV^ 












0.0095 GeV^ 








(OX0(3pll|)) 




0.107 GeV5 








(Oxo (3sM)) 




0.0022 GeV^ 









6. Conclusions 

As for the unpolarized J/4> yield, NRQCD factorization was consolidated at NLO 
by a global fit to the world's data of h adroproduction, photoproduction, two-photon 
scattering, and e"'"e~ annihilation,^^ which successfully pinned down the three CO 
LDMEs in compliance with the velocity scaling rules and impressively supported 
their universality. In a second step, NLO NRQCD predictions of J/i/j polariza- 
tion observables in various reference frames were confronted with measurements 
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in photoproduction at HERA and hadroproduction at the Tevatron and the LHC. 
In the case of hadroproduction at the Tevatron, the prediction of strongly trans- 
verse J/ip polarization in the helicity frame stands in severe contrast to the precise 
CDF II measurement which found the J/V' mesons to be unpolarized. Using 
the COLDME sets recently extracted from hadroproduction data by two other 
groupJ^^El 

does not help us to reach a satisfactory description of all the avail- 
able precision data. Thus, we conclude that the universality of the J/ip production 
LDMEs is challenged. Possible remedies include the following: 

(i) The eagerly awaited J/ip polarization measurements at the LHC might not 
confirm the CDF II results. 

(ii) Although unlikely, measurements at a future ep collider, such as the LHeC,S^ 
might reveal that the px distribution of J/ip photoproduction exhibits a drasti- 
cally weaker slope beyond px = 10 GeV, the reach of HERA, so that the LDME 
sets of Refs. \52\ [53l might yield better agreement with the data there. 

(iii) The assumption that the v expansion is convergent might not be valid for 
charmonium, leaving the possibility that the LDME universality is intact. 
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Fig. 1. NLO NRQCD filffSl compared to RHICpll TevatronpSEMZl LHCpSMEQEHSl] 
HERAPIMIMIlEP Hp], and KEK^ data. 



February 19, 2013 3:21 WSPC/INSTRUCTION FILE mpla 



12 Mathias Butenschoen, Bernd A. Kniehl 




PrlGeVj z E^[GeVl 



Fig. 2. ATLAS data on J/i/i inclusive hadroproductiorBSl and FTPS data on J/tp inclusive pho- 
toproduction in the fixed-target mode^l' compared to NLO NRQCD predictions evaluated using 
set A of CO LDMEs from Table [T] 
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Fig. 3. ZEUS data on J/xp inclusive photoproductiorlS] compared to NLO NRQCD predictions 
evaluated using set A of CO LDMEs from Table [T] 
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Fig. 4. The polarization param eter s A and u for direct photoproduction at HERA evaluated at 
NLO in the CSM and in NRQCeHSI using set B of CO LDME from Tableware compared to HiSSl 
and ZEud33] data. The theoretical uncertainties are due to scale variations in the CSM (blue 
bands) and include also the errors on the CO LDMEs (yellow bands) in NRQCD. 
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Fig. 5. The polarization parameters and for hadroproduction evaluated at NLO in the 
GSM and in NRQCC^QjSB using set B of CO LDME from Tableware compared to CDF data 
from Tevatron runs and and to ALICE dataH3 The theoretical uncertainties are due to 
scale variations in the CSM (blue bands) and include also the errors on the CO LDMEs (yellow 
bands) in NRQCD. 
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Fig. 6. The unpolarized J/t/; yields measured in e+e annihilation by BellefiSlin photoproduction 
by HIM and in hadroproduction by CDF itfZJ and ATL A&^Q as well as the J/tf> polarization 
observable Xg in the helicity frame as measured by CDF li^Slare compared with the NLO NRQCD 
predictions evaluated using the CO LDME sets of Refs.[T8l[52l|53] listed in Table[2] The theoretical 
errors in graphs a— g refer to scale variations, and those in graph d are obtained by also adding 
in quadrature the fit errors on the CO LDMEs according to Table [T] Graph h is taken over from 
Fig. 4 of Ref.[53l In graphs i-1, the central lines refer to the default CO LDME set of Ref.[52l and 
the theoretical errors are evaluated using the alternative CO LDME sets of Ref. 1521 



